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ABSTRACT 

We show that the nonlinear evolution of the cosmic gravitational clustering 
is approximately spatial local in the x-k (position-scale) phase space if the 
initial perturbations are Gaussian. That is, if viewing the mass field with modes 
in the phase space, the nonlinear evolution will cause strong coupling among 
modes with different scale k, but at the same spatial area x, while the modes at 
different area x remain uncorrelated, or very weakly correlated. We first study 
the quasi-local clustering behavior with the halo model, and demonstrate that 
the quasi-local evolution in the phase space is essentially due to the self-similar 
and hierarchical features of the cosmic gravitational clustering. The scaling of 
mass density profile of halos insures that the coupling between (x — k) modes 
at different physical positions is substantially suppressed. Using high resolution 
N-body simulation samples in the LCDM model, we justify the quasi-locality 
with the correlation function between the DWT (discrete wavelet transform) 
variables of the cosmic mass field. Although the mass field underwent a highly 
non-linear evolution, and the DWT variables display significantly non-Gaussian 
features, there are almost no correlations among the DWT variables at different 
spatial positions. Possible applications of the quasi-locality have been discussed. 

Subject headings: cosmology: theory - large-scale structure of the universe 
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1. 



Introduction 



The large scale structure of the universe was arisen from initial fluctuations through 
the nonlinear evolution of gravitational instability. Gravitational interaction is of long 
range, and therefore, the evolution of cosmic clustering is not localized in physical space. 
The typical processes of cosmic clustering, such as collapsing and falling into potential 
wells, the Fourier mode-mode coupling and the merging of pre-virialized dark halos, are 
generally non-local. These processes lead to a correlation between the density perturbations 
at different positions, even if the perturbations at that positions initially are statistically 
uncorrelated. For instance, in the Zel'dovich approximation (Zel'dovich 1970), the density 
field p(x, t) at (Eulerian) comoving position x and time t is determined by the initial 
perturbation at (Lagrangian) comoving position, q, plus a displacement S: 



The displacement S(q, t) represents the effect of density perturbations on the trajectories 
of self-gravitating particles. The intersection of particle trajectories leads to a correlation 
between mass fields at different spatial positions. Thus, the gravitational clustering is 
non-local even in weakly non-linear regime. 

On the other hand, spatial locality has been employed in the Gaussianization technique 
for recovery of the primordial power spectrum (Narayanan & Weinberg 1998). Underlying 
this algorithm is to assume that the relation between the evolved mass field and the initial 
density distribution is local, i.e. the high(low) initial density pixels will be mapped into 
high(low) density pixels of the evolved field (Narayanan & Weinberg 1998). Obviously, the 
localized mapping is difficult in reconciling the initially Gaussian field with the coherent 
non-linear structures, such as halos with scaling behavior. It has been argued that the 
locality assumption may be a poor approximation to the actual dynamics because of the 
non- locality of gravitational evolution (Monaco & Efstathiou 2000). Nevertheless, the 
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localized mapping is found to work well for reconstructing the initial mass field and power 
spectrum from transmitted flux of the Lja absorption in QSO spectra (Croft et al. 1998.) 
These results, joint with the data of WMAP, have been used to determine the cosmological 
parameters (Spergel et al 2003). However, the dynamical origin of the locality assumption 
remains a problem. It is still unclear under which condition the localized mapping is a good 
approximation. 

This problem has been studied in weakly non-linear regime under the Zel'dovich 
approximation. The result showed that the cosmic gravitational clustering evolution is 
spatially quasi-localized in phase (x-k) space made by the DWT decomposition (Pando, 
Feng & Fang 2001). In this approach, each perturbation mode corresponds to a cell 
in the x-k space, (x to x + Ax, k to k + Ak) with AxAk = 2n, and the density 
perturbation of the mode (x, k) is S(k, x). They demonstrated that, in the Zel'dovich 
approximation, if the initial density perturbations in each cells are statistically uncorrelated, 
i.e. {So(ki, xi)So(k2, x 2 )) oc k2 5^ 1X2 , where 5 K denotes for Kronecker delta function, 
the evolved 5(ki,x±) and S(k 2 ,x 2 ) will keep approximately spatially uncorrelated always, 
(5(ki, Xi)S(k 2 , x 2 )) oc 5^ , which is just what we call spatial quasi-locality in the dynamics 
of cosmic clustering. The spatial quasi-locality implies a significant local mode-mode 
coupling (different scales k at the same position x), but very weak non-local coupling 
between the modes (different positions x). The non-linear dynamical evolution is well 
developed along the direction of k-axis, rather than x-axis in the phase space. This 
quasi-locality has been justified in weak nonlinear samples such as the transmitted flux of 
QSO Lja absorption spectrum (Pando, Feng & Fang 2001). It places the dynamical base 
of recovering initial power spectrum from corresponding weakly evolved field via a localized 
mapping in phase space (Feng & Fang 2000). 

This paper is to extend the concept of quasi-locality to fully nonlinear regime. We 
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try to show that the quasi-locality of the cosmic clustering in phase space holds not only 
in the weak nonlinear regime, but also in nonlinear evolution. Since the nonlinear cosmic 
density field can be expressed by the semianalytical halo model (e.g. Cooray & Sheth 2002, 
and references therein), our primary interest is to study whether the quasi-locality could 
be incorporated in the halo model. We will first analytically derive the quasi-locality from 
the halo model, and then make a numerical test using high resolution iV-body simulation 
samples. 

The outline of this paper is as follows. §2 presents the statistical criterion of the 
quasi-local evolution of a density clustering in the x-k phase space. §3 shows that the 
density field evolution might be spatially quasi-localized in the phase space if the cosmic 
density field can be described by the halo model. Numerical tests on these predictions with 
N-body simulation samples are made in §4. Finally, the conclusions and discussions will be 
given in §5. 

2. Quasi-Locality in x-k Space 

2.1. DWT Variables of the Mass Field 

In physical space (x), the mode is Dirac delta function 8 D (x), and cosmic mass density 
field variable is p(x), while in scale space (k), the mode is the Fourier bases e lk x , and the 
field variable is p(k), which is the Fourier transform of p(x). In hybrid x-k phase space, one 
can use the complete and orthogonal bases of the discrete wavelet transform (DWT) as the 
mode function. The mass density field is then described by the DWT variables. 

Without loss of generality, we introduce the DWT variables by considering a density 
field p(x) in a cubic box of < x l < L, % — 1, 2, 3 and volume V = L 3 . We first divide the 
box into cells with volume L 3 /2 J1+ - 72+J3 , where j±, j 2 , js = 0,1.... For a given j = j 2 , j'3), 
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there are 2^+^ 3 cells labelled by 1 = (h,l 2 ,k), and k = 0, 1...2 ji - 1. The cell (h,l 2 ,k) 
occupies the spatial range lL/2 ji < Xj < (I + l)L/2 ji , % — 1,2, 3. Accordingly, indexes j and 
1 denote for, respectively, the scale and the position of the cells. In each dimension, we have 
Axi = L/2 ji and Aki = 2n/(L/2 j ), i.e. AxiAk = 2tt, or the volume of all cells in the x-k 
space is (2rr) 3 . 

Each cell (j,l) supports two compact functions: the scaling function 0j,i(x) and the 
wavelets Vj,i( x ) (Daubechies 1992, Fang & Thews, 1998, Fang & Feng, 2000). Both ^(x) 
and are localized in cell (j,l). The scaling functions 0j,i(x) are orthogonormal with 

respect to index 1 as 

/0j,i(x)0j,i'(x)dx = 5i,i'. (2) 

The scaling function 0j,i(x) is a low pass filter at cell (j,l). The scaling function coefficient 
(SFC) of the density field is defined by 

e J,i = / P(x)<Ai,i(x)dx, (3) 
which is proportional to the mean density in cell (j, 1). 

The wavelets ipj,\{x) are orthogonormal with respect to both indexes j and 1 

/ ^j,i( x M',i'( x ) rfx = ^j,Ai'- ( 4 ) 

The wavelets {?/>j,i(x)} form a complete and orthogonal base (mode) in the phase space. 
Therefore, the density field can be described by the wavelet function coefficients (WFCs) 
defined as 

ej,! = J dxp(x)V>j,i(x). (5) 

The WFCs e^i are the DWT variables of the density field. The DWT variables ej,i is the 
fluctuation of the density field around scales k = 27rn/L, n = (2 jl ,2 j2 ,2 js ) located at the cell 
1 with size Ax = (L/2- 71 , L/2^ 2 , L/2 J3 ). Since the wavelet ^j,i(x) is a band pass filter, in each 
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dimension, eji is a superposition of fluctuations filtered in the waveband k ± Ak/2, where 
k = 2ir2 j /L, and AA; = 2ir/ Ax = k. This decomposition of the fluctuation is optimized in 
sense that the size of cell Ax is adaptively chosen to match with the perturbations at a 
given wavenumber. 

The DWT bases generally have vanishing moments, i.e. 

J xVj,i(x)dx=0, (6) 

where t — 0, 1...M — 1 and M > 1. M is dependent on wavelet. For the wavelet Daubechies 
2n, we have M = n. Thus, the Fourier transform of wavelet V^k) = / V'j.iW exp(— ik-x)rfx 
has a compact support in the wavenumber space {k}. From eq.(6) we have / /w/>j,i(x)<ix = 
for all j, 1, and thus, 

ej,i = pj dx<5(xM,i(x) (7) 

where the density contrast <5(x) = [p(x) — p]/p, p is the mean density. Since (5(x)) = 0, we 
have also 

<q,i> = o (8) 

Because the set of wavelets is complete, give a complete description of the density 
field p(x), i.e. one can reconstruct p(x) or 5(x) in terms of variables e^i as 

2.2. Quasi-Locality of Gaussian Fields 

The initial density perturbation of the universe <5(x, tj) is believed to be a Gaussian 
random field with correlation matrix of the Fourier variables <5(k, ti) 

(5(k,t l )5(k\t l )}=P(k,U)5«, (10) 
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and all higher order cumulant moments of <5(k, ti) vanish. The function P(k, tj) of eq.(10) 
is the initial power spectrum. The Kronecker delta function 5^ k , in eq.(10) indicates that 
the initial perturbation for each mode k is independent, or localized in /c-space. 

Generally, if the initial Fourier power spectrum, P(k, tj), is colored, i.e. A;-dependent, 
the correlation matrix of variables other than the Fourier mode will no longer be diagonal. 
For instance, correlation function in x-space will be (<f(x,i i )<J(x',* i )) = £(|x-xr,ti), 
which is the Fourier counterpart of P(k, ti). However, the correlation matrix of the DWT 
variables of a Gaussian field is always diagonal or quasi-diagonal, regardless the Fourier 
power spectrum P(k) is white or colored. That is, the correlation function of e^i is localized 
with respect to (j, 1) as 

%AU)iyAU)) - PiAW&x,, (11) 
where in eq.(ll) is the DWT power spectrum of the field. 

The reason of the diagonality of eq.(ll) is as follows. First, the WFC ej 5 i is given by 
a linear superposition of the Fourier modes 5(k) in the waveband around fcj ~ 2n2 ji /L 
(i = 1,2,3), with different j consist of the Fourier modes <5(k) in different k bands. 
While for Gaussian fields, the Fourier modes in different wave bands are uncorrelated in 
general [eq.(10)], and therefore, there might be no correlation between the DWT modes of 
j and j' if j 7^ j'. This yields the quasi-locality of <5jj,. Second, the phases of the Fourier 
modes of Gaussian field are independent and random. For a superposition of the random 
phased Fourier modes <5(k) in the band from k to k + Ak, the spatial correlation length can 
not be larger than that given by the uncertainty relation Ax ~ 2ir/ Ak ~ 2~ j L. Moreover, 
the non-zero regions of two DWT modes ip^i and ip^r with (1 ^ 1') have spatial distance 
A|xj| ~ 2~i i L\l i — Consequently, all off-diagonal elements (1 ^ V) vanish or are much 
smaller than diagonal elements, i.e., 

\Jdxf dxVj,i(x)^(x - x', ti)il>y iV (x.') 



Jdx/dxVj,i(x)e(x-x',^)Vj,i(xO 



«1, if 1^1' or j^j'. (12) 
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Thus, the correlation function of the DWT variables of a Gaussian field is rapidly decaying 
when |1 — l'| > 1 and |j — j'| > 1. The diagonal correlation function described by eq.(ll) is 
a generic feature of a Gaussian field in the DWT representation. 

It should be pointed out that the WFC correlation function eq.(ll) is different from 
the ordinary two point correlation function (5(x)5(x')) = £(x, x'). The former is the 
correlation between two modes in (x — k) phase space, while the later is for two modes in 
x-spa.ce. Explicitly, eq.(ll) describes the correlation of perturbation modes in the waveband 
k — > k + Ak between positions 1 and 1', and so, it is sensitive to the phases of modes. The 
ordinary two point correlation is not sensitive to the phase of perturbations. In the DWT 
analysis, an analogue of the ordinary two point correlation function is defined by by the 
correlation between SFCs, i.e. (cj,\(ti)ey,\'(ti)) . Since the scaling function (j>- h \ is a low-pass 
filter on scale j and at position 1, the correlation function {e^i(ti)eyy(ti)) behaves in a 
similar way as (<5r(x)<5r(x')), where <5r(x) is a filtered density field smoothed on the scale 
R ~ 2~ j L. However, as ^j.i is a high-pass filter, the WFC covariance (ej,iej',i') shows quite 
different statistical features from the SFC correlation, e.g., it is always quasi-diagonal or 
even fully diagonal for a Gaussian field. Generally, it has been shown for many analytically 
calculable random fields that the SFC correlation is significantly off-diagonal, while the 
WFC correlation is exactly diagonal (Greiner, Lip & Carruthers, 1995) 

Moreover, within a given volume in the x-k space, such as Vd 3 k, the number of the 
DWT modes {j,l} is the same as that of the Fourier modes k. Accordingly, P(k, tj) can 
be expressed as a linear superposition of Pyi(ti), and vice versa. Equivalently, the Fourier 
power spectrum can be replaced by the DWT power spectrum (Fang & Feng, 2000). 
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2.3. Statistical Criterions of Quasi-Locality 

If the evolution of the comic mass field is localized, the evolved density field <5(x) 
at a given spatial point is determined only by the initial density distribution <5(x, tj) at 
the same point. As emphasized in §1, this locality is inconsistent with the non-local 
behavior of gravitational clustering. However, the evolution of the comic mass field can be 
quasi-localized in sense that the correlation between the DWT variables of the evolved field 
is always spatial diagonal if the initial correlation function is diagonal, such as eq.(ll). For 
perturbation modes in a waveband k — > k + Ak, the quasi-localized range is Ax ~ 2n/ Ak. 

A quasi-localized evolution means that the auto-correlation function of the DWT 
variables is always diagonal, or quasi-diagonal (£j,i£y,\) when 1 ^ 1', if it is 

diagonal initially. Thus, one may place a statistical criterion for the quasi-locality as 

«jj,(A1) « 1, ifAl^O, (13) 

where 

^' (A1) = <wn\hW- (14) 

This is a normalized correlation function of the DWT modes, i.e. ftjj(Al = 0) = 1. 

The auto-correlation function of the DWT variables, measures the correlations 

between the perturbation modes on scales j = (ji, j 2 , jz) and j' = j^fz) at two cells 
with a vector distance given by (hL/2 jl - l' 1 L/2 j 'i,l 2 L/2 j2 - l' 2 L/2 j 2,l 3 L/2 js - l' 3 L/2 j *). In 
the case of Al = 0, ftjj'(O) gives the correlation between fluctuations on scale j and j' at 
the same physical area. Therefore, if condition (13) holds for all redshifts, the dynamical 
evolution of the mass field is basically spatial localized in the DWT bases. Comparing the 
condition eq.(13) with eq.(12), we see that the cosmic field undergoing a local evolution is 
different from its Gaussian predecessor by the factor 6-^,, but not 5^ v . In other words, the 
evolution leads to the significant scale-scale (j, j') coupling, rather than modes at different 
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locations 1 7^ 1'. 

One can also construct the criterions for the quasi-locality using higher order 
correlations among the DWT variables e^. For instance, a (p + q) order statistical criterion 
is given by 

CJ3?(A1 # 0) ~ 1. (15) 

where 

= w&' (16) 

where p and q can be any even number. Obviously, C?y ~ 1 for Gaussian fields. C?y(0) 7^ 1 
corresponds to a local scale-scale correlation, while C?y(A\ ^ 0) ^ 1 a nonlocal scale-scale 
correlation. The quasi-local evolution of cosmic mass field requires that nonlocal scale-scale 
correlation is always small. 

It should be pointed out that the statistical conditions eqs.(13)-(16) are not trivial 
because the DWT basis does not subject to the central limit theorem. If a basis subjects 
to the central limit theorem, the corresponding variables will be Gaussian even when the 
random field is highly non-Gaussian. In this case, eqs.(13) and (15) may be easily satisfied, 
but it does not imply that the evolution is localized, or quasi-localized. Statistical measure 
subjected to the central limit theorem is unable to capture non-Gaussian features of the 
evolution. 

3. Quasi-Local Evolution in Halo Model 

We will show, in this section, that the statistical criterions of §2.3 are fulfilled if the 
cosmic mass field can be described by the halo model. 
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3.1. 



The Halo Model 



The cosmic clustering is self-similar and hierarchical, as the dynamical equations of 
collisionless particles (dark matter) do not have preferred scales, and admits a self-similar 
solution as well as the initial density perturbations are Gaussian and scale free. The halo 
model further assumes that all mass in a fully developed cosmic mass field is bound in halos 
on various scales (Neyman & Scott 1952, Scherrer & Bertschinger 1991). Thus, the cosmic 
mass field in non-linear regime is given by a superposition of the halos 



p(x) = YlPii* - x i) = Xl m * M ( x ~ x i> m *) = / rfx'^mi5 D (x' - Xi)ti(x - x',TOi), (17) 

i i i 

where p^(x — Xj) is the density profiles of halo with mass at position Xi, and u(x — x;, m,) 
is the density profile normalized by / rfxpj(x — Xj) = / dxrajw(x — Xj, m,) = mj. 

There are several different versions of the halo density profiles, such as 
m(x, ITLi ) oc l/(r/r;) a [l + (r/r i )] f3 with a = 1,(3 = 2 (Navarro, Frenk & White 
1996), and u(x,mj) oc l/(r/rj) a [l + (r/^)^] with a = 3/2,(3 = 3/2 (Moore et al. 1999). A 
common feature of the halo density profiles is self-similar, which implies that the indexes 
a and (3 should be mass-independent. The mass dependence is only given by r^, which 
characterizes the size of the mj halo. The details of the profiles are indifferent for the 
problem we try to study below. What is important for us is only that one can set a 
self-similar upper limit to the normalized halo density profile as 



where r = |x|, C is a constant and the index 7 is mass-independent. The m^-dependence of 
Ti are not stronger than a power law as r% oc raj. 

The halo model also assumes that the halo-halo correlation function on scales larger 
than the size of halos is given by the two-point correlation functions of the linear Gaussian 




(18) 
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field with a linear bias correction, or by the correlation functions of quasi-linear field. 
Therefore, no assumption about higher order halo-halo correlations on large scales is needed. 

In this model, the time-dependence of the field is mainly given by the mass function 
of halos, n(m,t), which is the number density of the halos with mass to at time t. In 
the Press & Schechter formalism (1974), the mass function is determined by the power 
spectrum of the initial Gaussian density perturbation. Moreover, the self-similarity of the 
halo density profiles insure that eq.(18) holds for all time. The cosmic evolution only leads 
to the parameters on the r.h.s of eq.(18) to be time- dependent. 



3.2. Quasi-Locality of the DWT Correlation Function 

With eq.(17), the DWT variable of the cosmic mass field in the halo model is given by 
e j5 i = J dxp(x)V>j,i(x) = J2 m i J rfx «( x - x i ,m i )^j i i(x). (19) 
The auto-correlation function of the DWT variables is then 

(MM') = %M,r) h + (Wy,r) h \ (20) 

where the first and second terms on the r.h.s. are usually called, respectively, 1- and 2-halo 
terms. They can be written in the explicit form, 

(£j,i£y,\') h = J dmn{m)m 2 J c?x 1 -u(x 1 ,m) (21) 
J dx J dx / '0j ! i(x)u(x 1 + x — x', m^j/^x') 

(£},\£y ,\') hh = J dmin(mi)mi J rfm 2 n(m 2 )m 2 J rfxi J rfx 2 -u(xi, mi)u(x 2 , m 2 ) (22) 
J dx J rfxVj,i(x)^(x-x' + Xi - x 2 ,mi,m 2 )V' j / ) i/(x / ), 

where n{m) = (^2i5 D (m — mj)<5 D (x — x;)) is the number density of halos with mass to, and 
£(x - x', toi, to 2 ) = (J2i^j 5 D (mi - rrii)5 D (x. - x ; )5 d (to 2 - m i )5 D (x l - x ; )) / 'n(mi)n(m 2 ) is 
the halo-halo correlation function. 
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We show now that the DWT correlation function (ej,iej',i') is quasi-diagonal, or fast 
decaying with respect to the spatial distance 1 — 1'. First, consider the 2-halo terms eq. 
(22). According to the halo model, the two-point correlation function £(xi — x 2 ,mi,m 2 ) is 
determined by the linearly Gaussian density field. In fact, as having been discussed in §2.2, 
the DWT correlation function of a Gaussian field is generally diagonal [eq.(12)], regardless 
the Fourier power spectrum is white or colored. Therefore, the DWT correlation function 
(£j,iej',i') contributed from the 2-halo term should be quasi-local. 

The DWT integral in the 2-halo term of eqs.(22) is not completely the same as eq.(12), 
as the halo-halo correlation function in eq.(22) is £(x — x' + Xi — x 2 ), while in eq.(12) 
is £(x — x'). For halos with size less than the scale j considered, the factor |xi — x 2 | in 
correlation function £ hh is smaller than the size of the DWT mode, L/2- 7 , and so, the factor 
X! — x 2 can be ignored in comparison with the variable x — x'. Moreover, by definition 
of the halo model, the correlation functions £ hh does not include the contributions from 
halos with sizes larger than |x — x'|. Thus, the cases of |xi — x 2 | > |x — x'| can always be 
ignored. The 2-halo term essentially follows eq.(12), and is always approximately diagonal 
with respect to spatial indexes 1 and 1'. 

To analyze the 1-halo term eq.(21), we use the following theorem of wavelets (Tewfik 
& Kim, 1992). Because the DWT bases is self-similar, for any 1-D power law function 

f(x) oc x -7 , we have 



where ipj,i{x) is wavelet in 1-D space, and C is a constant. Therefore, while using wavelets 
with large enough M [eq.(5)], the integral eq.(23) is quickly decaying with the spatial 
distance \l — In other words, besides two nearest position of I — V — ±1, there is no 
correlation between modes at different position I ^ V . 



The 3-D integral J(|l-l'|) = | / dxj dx / Vj,i( x ) u ( x i + x - x', m)^j,i'( x ')l in eq.(21) 




(23) 
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has the similar structure as eq.(23). Using the upper limit eq.(18), one can estimate the 
decaying of the integral — l'|) with |1 — l'|. When |x x | < |x — x'|, we can expand the 
function |xi + x — x'|~ 7 in terms of |xi|/|x — x'|. Thus, applying eq.(23) term by term, 
we have the decaying behavior of I at least as fast as 2~ jl+j2+j3 \\ — l'|-7- 2M . Similarly, 
in the case of jx^ > |x — x'|, one can expand the function |x x + x — x'|~ 7 by the factor 
|x — x'|/|xi|. Following from eq.(23), the decaying behavior of / is then decreasing with 
|1 — l'| as 2 _ - 7l+ - 72+ - 73 |l — l'| _4M . Accordingly, the 1-halo term in the correlation between 
modes j, 1 and j, 1' will generally decay as |1 — 1'|-t- 2M or |1 — l'|~ 4M . The correlation 
function (ej,iej,i')' i is then approximately diagonal w.r.t. the spatial index 1 and 1'. This 
result is largely valid due to varying with m 8 by a power law. 

Proceeding in the similar way as above, we can also show the diagonality of the 
correlation function %,i€y,r) h - In this case, instead of eq.(23), we use the following theorem 



Consider that the cell (j, /) has the same physical position as cell (j', I') if V — I x 2 j '~ j , 
the theorem eq.(24) also yields that the 1-halo term of the two modes j,l and j',1' will 
decline with the physical distance between the two modes as |1 — (Y /2 j '~ j )\~' Y ~ 2M or 
|1_ (y/2i'-i)\-*M . 

Based on above discussions, one may draw the conclusion that if the halo model 
is a good approximation to the cosmic density field p(x) in the nonlinear regime at all 
time, their correlation matrix in the DWT representation remains quasi-diagonal forever, 
and the evolution is quasi-local. This result is based on the self-similarity of the density 
profiles of halos and the weakly nonlinear correlation between halos. The self-similar scaling 
ensures that the non-local correlation among the DWT variable is uniformly suppressed, 
independent of the mass of halos. Mathematically, the correlations between the DWT 



(Tewfik & Kim, 1992) 



/ 
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if f > j- 



(24) 
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modes of perturbations at different physical places are uniformly converging to zero with 
the increasing of M if the index 7 is mass-independent. 



3.3. Quasi-Locality of Higher Order Statistics 

To show the quasi-locality of higher order statistical criterion, we use the hierarchical 
clustering or linked-pair relation (Peebles 1980), which is found to be consistent with the 
halo model. For the 3 rd order correlations, the linked-pair relation is 

(^."^^(^i-^ia) * g 3 [(5(x 1 )5(x 2 ))(5(x 1 )5(x 3 )) (25) 

+ 2 terms with eye. permutations], 

where the coefficient Q 3 might be scale-dependent. Subjecting eq.(25) to a DWT by 3 rd 
order basis ^j 1 ,i 1 (xi)^j 2 ,i 2 (x 2 )^j3,i 3 ( x 3) ) we have 

(^1,11^2,12^3,13) - Q3 XI a ji,iij',i'a",i"[( ? j',i' ? j2,i2)( ? j",i"ej3,i 3 ) ( 26 ) 
+ 2 terms with eye. permutations]. 

where aj\ {/1 ^ is given by the 3-wavelet integral, 

4,ii;jM';j",i" =/^ji,ii( x )^jM'( x )^j'M"( x )* c - ( 27 ) 

Since V'j.iM is localized in the cell (j,l), a 3 is significant only if the three cells (ji,li), 
(j', 1') and (j", l") coincide with each other at the same physical area. Thus, by virtue of the 
locality of correlations (cj/ i'ej 2 ,i 2 ) and i"ej 3 ,i 3 ) (§3.2), it is easy to see that (ej li i 1 ej 2) i 2 ej3 i i 3 ) 
is small if the cells (ji,li), (j2, 12) and ^3,13) are disjoint in the physical space. Since Q 3 
does not depend on x, the result of locality will keep valid when Q 3 is scale-dependent, 

Obviously, the 3 nd order result can be generalized to n th order DWT correlation 
function. The integral of n wavelets V!ji,ii> V'j2,i2v'0j„,in i s zero or ver y small otherwise the 
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n cells (ji,li), (j2 5 12) •• • (jn; l n ) coincide in the same physical area. In addition, all terms 
on the r.h.s. of the hierarchical clustering relation consist of linked 2 nd DWT correlation 
function, the 11 th order DWT correlation function (ej 1 ,i 1 ej 2) i 2 ...ej ni i n ) should be localized. 
Thus, the criterion eq.(15) and other higher order criterions will be satisfied in general. 

In summary, if the cosmic density field is evolved self-similarly from an initially 
Gaussian field, the spatial quasi-locality is true at all time, i.e. 

• the second and higher order correlation functions of the DWT variables, ej,i of the 
evolved field is quasi-diagonal with respect to the position index 1. 

For those type of fields, the possible non-Gaussian features with the DWT variables are 
mainly 

• non-Gaussian one-point distribution of the DWT variables e^; 

• local scale-scale correlation among the DWT variables. 

Above three points are the major theoretical results of this paper. 

4. Testing with N-body Simulation Samples 
4.1. Samples 

To demonstrate the quasi-locality of the evolved cosmic density field, we use samples 
produced by Jing & Suto (2002). The samples are given by high resolution N-body 
simulation, running with the vectorized-parallelized P 3 M code. The cosmological model was 
taken to be LCDM model specified by parameters (Q ,A,a 8 ,r) = (0.3,0.7,0.9,0.2). The 
primordial density fluctuation is assumed to obey the Gaussian statistics (this is important 
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for us), and the power spectrum is of the Harrison-Zel'dovich type. The linear transfer 
function for the dark matter power spectrum is taken from Bardeen et al (1986). 

The simulation was performed in a periodic, cubical box of size 100h~ 1 Mpc with a 
1200 3 grid points for the Particle-Mesh (PM) force computation, and 512 3 particles. The 
short-range force is compensated for the PM force calculation at the separation less than 
e = 2.7H, where H is the mesh cell size. The simulations are evolved by 1200 time steps 
from the initial redshift = 72. The force resolution is 20 h" 1 kpc for the linear density 
softening form. It is noted that our statistic tests are performed on scales > 0.2 h" 1 Mpc. 

We have one realization. In the practical computations, we divide the 100 h -1 Mpc 
simulation box into 8 subboxes each with size L = 50h _1 Mpc. Accordingly, the ensemble 
average and la variance are obtained from those 8 subboxes. 

4.2. Two-Point Correlation Functions 

Before showing the quasi-locality, we first calculate the correlation function of the 
SFCs, i.e. (eyey/). As has been discussed in §2.2, the correlation function (eyey) is 
actually analogue to the ordinary two-point correlation function £(r), where r = x — x'. 
Since ey is a filtered density field smoothed by the scaling function on the scale j, it is 
expected that the correlation function (eyey) will display similar feature as £(r). 

Figure 1 presents the r-dependence of (eyey) at j = (7,7,7), corresponding to the 
smoothed field filtered on the linear scale 50/2 7 = 0.39 h _1 Mpc. The spatial distance 
between the cells 1 and 1' is r = |1 — l'|50/2 7 h _1 Mpc. In this calculation, we applied 
wavelet Daubechies 4 (Daubechies, 1992). As expected, the r- or |1 — 1' (-dependence of 
the correlation function (eyey/) shows the standard power law, (eyey/) oc r~ a with the 
index a ~ 2. That is, this correlation function is not localized. Meanwhile in Fig.l, the 
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contributions of the 1-halo and 2-halo terms are also plotted, respectively. The 1-halo 
term is calculated from eq.(21) with the NFW density profiles of halos, and the 2-halo 
term is given by eq.(22) in which the correlation function £(r) is given by the linear power 
spectrum. The nonlinear clustering is largely due to the 1-halo term. Therefore, one can 
see that the nonlinear evolution of gravitational clustering cause correlations on scales of 
several Mpc. 

4.3. Justifying the Quasi-Locality 

We now study the quasi-locality of the clustering with the correlation function (ej,iej',i'), 
which is used in the criterion Kjj/(A1) [eq.(13)]. First, we take the same parameter as Fig. 
1, j = j' = (7, 7, 7), and r = |Al|50/# = |1 - l'|50/2- ? ' h" 1 Mpc, and also we used wavelet 
Daubechies 4 (Daubechies, 1992), which has M — 2. The result of DWT correlation 
function (ej,iej,i')- is shown in Fig. 2. The solid circle at r = in Fig. 2 corresponds to 
Al = 0, or 1 = T, and other solid circles from left to right correspond, successively, to 
Al= 1, 2, 3... 

From Fig. 2, we can see immediately that the shape of the r- or | Al|-dependence of 
(£j,iej,i') is quite different from the "standard" power law. The correlation function (ej,i£j,i') 
is non-zero mainly at point r = or |A1| = 0. At |A1| = 1, the correlation function 
(ej,iej,i±i) drops to tiny values around ~ 0. For |A1| > 1, the correlation function basically is 
zero. The correlation length in terms of the position index 1 is approximately zero, namely, 
the covariance is diagonal. This is the spatial quasi-locality. In Fig. 2, we also 

plot the 1- and 2-halo terms, an d (^j,i^j,v) hh ■ Although 1-halo term is dominated 

by massive halos, the covariance (ej j iej i i/) h is also perfectly quasi-localized because of the 
self-similarity of density profiles of massive halos. The 2-halo term is zero at r = 0, because, 
by definition, eq.(22) does not contain the contribution of autocorrelations of halos. 
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Figure 3 presents «jj/(Al) vs. r, for j = j' = (4,4,4), (5,5,5) and (6,6,6). The 
physical distance r is the same as Fig. 2, given by r = |A1|50/2 J = |1 — l r 1 50 / 2 J h -1 Mpc. 
The solid circle at r = corresponds to |A1| = 0, at which, by definition, we have the 
normalization Kjj(O) = 1. Other points from small to large values of r correspond to Al = 
1,2,.... successively. Clearly, all /tjj(Al) for Al > are less than 10 -6 , which is actually 
from the noises of sample. The result implies that for all calculable points of r > 0, the 
correlation is negligible, and satisfies the criterion eq.(13). 

We also calculated k^i (Al) for modes of j = j' = (j±, j 2 , js), but j\ 7^ 3i 7^ 33- Most of 
these cases shows Kjj/(A1) ~ if |A1| > 0. Only exception is for the cases of ji, j 2 < js, 
and |A1| = |/ 3 — l' 3 \ = 1. As an example, Fig. 4 presents Kjj(Al) for modes j = (5,5,6), 
and r = |1 3 — 1 3 |50/2 J ' 3 . It shows that (ej,iej,i') = if |Z 3 — l' 3 \ > 1. The non-zero value at 
1^3 ~~ ^3 1 = 1 1S about 20% of that at Al = 0. This result is consistent with the theorem (23), 
which only requires the suppression for modes \l — l'\ > 1, but may not work for \l — l'\ = I. 
It implies that the clustering may give rise to the correlations between nearest neighbor 
cells in phase space. However, the cell resolved by ji,j 2 < J3 is a rectangle in the physical 
space, and the shortest edge is given by j 3 , the physics distance of |/ 3 — l' 3 \ = 1 is still less 
than whole size of the rectangle. Thus, it could be concluded that the covariance (ej,iej,i') 
is always quasi-diagonal in the sense that all members with 1 and 1' are almost zero if the 
distance r between 1 and 1' is larger than the size of the cell (j±, j 2 , j 3 ) considered. In Fig. 
4, we show also a result calculated with wavelets Daubechies 6 (D6), for which M — 3. It 
yields about the same results as D4. 

For the correlation between modes (j,l) and (j', 1') with j ^ j', we use the criterion 
of eq.(14). In this case, the physical distance between two cells is r = |r| and 
rj = |/j/2 Ji — l' i /2 : >'i\L h -1 Mpc. Figure 5 plots «jj/(Al) vs. r for modes j = (j,j,j), 
j' = (j + 1, j + 1, j + 1), and j = 3,4,5 and 6. All the values of «jj/(Al) in Fig. 5 are 
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not larger than 1CT 5 , and are much less than the diagonal terms (ef i), or (e? ;1 /)- We found 
this result is generally true for all the cases of j 7^ j'. That is, the second order correlation 
between two modes with different scales j and j' are always negligible, regardless the 
indices 1 and 1'. In other words, the covariance of the WFC variables (ej,iej',i') generally are 
quasi-diagonal. 

As for the high order statistics by criterion eqs.(15) and (16), we can cite some previous 
calculations of the non-local scale-scale correlation defined by 

Cff(Al) = ,^| +A1> , . (28) 

\ e j,lA e j,l+Al/ 

which is the criterion eq.(15) with p = q = 2. It has been shown that either for the APM 
bright galaxy catalog (Loveday et al. 1992) or mock samples of galaxy survey (Cole et 
al. 1998), the non-local scale-scale correlation always yields |Cj 2 j 2 (Al) — 1| <C 1 if Al > 
(Feng, Deng & Fang 2000). Although this work was not for addressing the problem of the 
quasi-locality, the result did support the quasi-locality up to the 4 th order statistics. 



4.4. Non-Gaussianity Revealed by DWT Variables 

As discussed in §2.3, it is necessary to show that the random variables e^i are 
non-Gaussian for evolved fields. The possible non-Gaussian features with the DWT 
variables are (1) non-Gaussian one-point distribution of ey, and (2) local scale-scale 
correlations (§3.3). 

In Fig. 6, we plot the one-point distribution on scales j = (j, j, j, ) and j = 5, 6, 7 and 
8. It illustrates that the kurtosis of the one-point distribution is high. The PDF (probability 
distribution function) is approximately lognormal. The 4 th order local scale-scale correlation 
Cg 2 (Al = 0) is plotted in Fig. 7. It shows Cj^Al = 0) > 1 on small scales (j = 5, 6 
and 7), while random data gives C?f (Al = 0) = 1 on all scales. The evolved field is highly 
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non-Gaussian, although it is always quasi-localized. 

5. Discussions and Conclusions 

We showed that the cosmic clustering behavior is quasi-localized. If the field is viewed 
by the DWT modes in phase space, the nonlinear evolution will give rise to the coupling 
between modes on different scales but in the same physical area, and the coupling between 
modes at different position is weak. The quasi- local evolution means that, if the initial 
perturbations in a waveband k ± Ak/2 and at different space range Ax is uncorrelated, the 
evolved perturbations in this waveband at different space range Ax will also be uncorrelated, 
or very weakly correlated. In this sense, the nonlinear evolution has memory of its initial 
spatial correlation in the phase space. This memory is essentially from the hierarchical 
and self-similar feature of the mass field evolution. The density profiles of massive halos 
obey the scaling law [eq.(18)], and therefore, the contributions to the non-local correlation 
function from various halos are uniformly suppressed. 

It has been realized about ten years ago that some random fields generated by a 
self-similar hierarchical process generally shows locality of their auto-correlation function in 
the phase space, if the initial field is local, like a Gaussian field (Ramanathan & Zeitouni, 
1991; Tewfik & Kim 1992; Flandrin, 1992). Later, this result are found to be correct 
for various models of structure formations via hierarchical cascade stochastic processes 
(Greiner et al. 1996, Greiner, Eggers & Lipa, 1998). These studies implies that the local 
evolution and initial perturbation memory seems to be generic of self-similar hierarchical 
fields, regardless the details of the hierarchical process. It has been pointed out that models 
for realizing the self-similar hierarchical evolution of cosmic mass field, such as the fractal 
hierarchy clustering model (Soneira & Peebles 1977), the block model (Cole & Kaiser 
1988), merging cell (Rodrigues & Thomas 1996), have the same mathematical structures 
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as hierarchical cascade stochastic models applied in other fields (Pando et al 1998, Feng, 
Pando & Fang 2001). Obviously, the local evolution can be straightforward obtained in 
those models. 

The DWT analysis is effective to reveal the quasi-locality in phase space. Such 
quasi- locality is hardly described by the Fourier modes 5(k), as the information of spatial 
positions is stored in the phases of all Fourier modes. Moreover, the Fourier amplitudes 
|<5(k)| subject to the central limit theorem, and are insensitive to non-Gaussianity. The 
wavelet basis, however does not subject to the central limit theorem (Pando & Fang 1998), 
which enable us to measure all the quasi-local features with the statistics of e^i. 

The quasi-locality of the DWT correlation is essential for recovery of the primordial 
power spectrum using a localized mapping in phase space. Such mapping has been 
developed in recovering the initial Gaussian power spectrum from evolved field in the 
quasi-linear regime (Feng & Fang 2000). By virtue of the quasi-locality in fully developed 
fields, we would be able to generalize the method of localized mapping in phase space to 
highly non-linear regime. 

The quasi-local evolution may also provide the dynamical base for the lognormal 
model (Bi, 1993; Bi & Davidsen 1997, Jones 1999). The basic assumption of the lognormal 
model is that the non-linear field can be approximately found from the corresponding linear 
Gaussian field by a local exponential mapping. The local mapping is supported by the 
quasi-local evolution. We see from Fig. 6 that the PDF of evolved field is about lognormal. 
Therefore, in the context of quasi-local evolution, a local (exponential) mapping from the 
linear Gaussian field to a lognormal field might be a reasonable sketch of the nonlinear 
evolution of the cosmic density field. 
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Fig. 1. — The SFC correlation function (ej,iej,i') vs. r for the simulation data (hexagon), 
1-halo (square) and 2-halo (circle) terms. The scale j is taken to be (7,7,7). The physical 
distance is given by r = |1 — l / |50/2- 7 h -1 Mpc. 
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Fig. 2. — The DWT variable (WFC) correlation function (ej,iej,i') vs. r for the simulation 
data (hexagon), 1-halo (square) and 2-halo (circle) terms. The scale j is taken to be (7, 7, 7). 
The physical distance is given by r = |1 — l / |50/2- ? h _1 Mpc. 
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Fig. 3. — /tjj vs. r for the simulation data, j are taken to be (4,4,4) (top panel), (5,5,5) 
(middle), and (6,6,6) (bottom). The physical distance is given by r = |1 — l'lSO/2- 7 h^ 1 Mpc. 
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Fig. 4. — Kjj'(AZ) vs. r for the simulation data. The scales are j = (5,5,6) The physical 
distance r = |Al|50/2 6 h" 1 Mpc. Wavelets D4(M = 2) and D6(M = 3) are used. 
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Fig. 5. — /tjj + i(Al) vs. r for the simulation data. The scale indexes are taken to be 
j = (j, 3, 3, ), i' = (j + 1, 3 + 1, 3 + 1), and j = 3, 4, 5 and 6, where Al = 1 - V/2. The physical 
distance r = |A1|50/2 J ' h" 1 Mpc. 
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Fig. 6. — The one-point distribution of ej,i for j = ) and j = 5, 6, 7 and 8. 
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Fig. 7. — The 4 th local scale-scale correlation Cj 2 j +1 (A1 = 0) vs. j for the simulation data 
(solid circle) and Gaussian random sample (circle). 



